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Extensions of the standard model (SM) with gauge symmetry SU{'i)c ® SU{A)l ® U{l)x (3- 
4-1 extensions) where anomaly cancellation takes place between the fermion families (three-family 
models) predict the existence of two new heavy neutral gauge bosons which transmit flavor changing 
neutral currents (FCNC) at tree-level. In this work, in the context of a three-family 3-4-1 extension 
which does not contain particles with exotic electric charges, we study the constraints coming from 
neutral meson mixing on the parameters of the extension associated to tree level FCNC effects. 
Taking into account experimental measurements of observables related to K and B meson mixing 
and including new CP-violating phases, we study the resulting bounds for angles and phases in 
the mixing matrix for the down-quark sector, as well as the implications of these bounds for the 
' modifications in the amplitudes of the clean rare decays if^ ~> -k^vv, Kl — > ■k'^vv, Kl tv''^1^1~ 

= e,/i) and Btj/s -5> ^"""M"- 

> ' 

O . PACS numbers: 12.60.Cn, 13.20.Eb, 13.20.He 

^ I. INTRODUCTION 

' ' \ Flavor violating couplings of ordinary fermions to extra neutral gauge bosons and to new scalar fields arise in many 
extensions of the standard model (SM). Two simple and interesting examples are the 3-3-1 model in which the SM 
I ■ gauge symmetry is enlarged to SU{3)c SU{3)l (8) U(\)x and the 3-4-1 extension where the enlargement is 

done to the gauge group SU (3) c^E) SU (4) l ^ U{l)x [SHlOj. These extensions share the important feature of addressing 
f-H ' the problem of the number of fermion families in nature in the sense that anomaly cancellation among the families 
I— I, (three- family models) takes place only if we have an equal number of left-handed triplets and antitriplets (in the 
3-3-1 model) or an equal number of 4-plets and 4*-plets (in the 3-4-1 extension), taking into account the color degree 
of freedom. As a consequence, the number of fermion families Nf must be divisible by the number of colors Nc of 
• SU{3)c- Moreover, since SU{3)c asymptotic freedom requires Nc < 5, it follows that Nj ^ Nc = 3. Cancellation of 
' chiral anomalies among the families thus forces one family of quarks, in the weak basis, to transform differently from 
the other two under SU{4)l <E) U{1)x- As a result, the chiral couplings of quarks to the new neutral gauge bosons 
' are, in general, nonuniversal and, when rotating to the quark mass eigenstate basis, flavor changing neutral currents 
' (FCNC) induced by fermion mixing arise. 

In particular, the 3-4-1 extension predicts the existence of two heavy neutral gauge bosons Z' and Z" which, in 
general, mix up with the ordinary Z boson of the SM. In contrast with the SM where FCNC processes occur only at 
the loop level, in the 3-4-1 extension these new gauge bosons can transmit FCNC at tree level and, consequently, the 
study of these processes can set stringent bounds on the Z' and Z" masses and mixing. Moreover, since in general 
each flavor couples to more than one Higgs 4-plet, FCNC coming from the scalar sector can also be present. 

For simplicity, in this work we will restrict ourselves to 3-4-1 extensions without exotic electric charges in the 
fermion sector, that is, without electric charges different from ±2/3 and ±1/3 for exotic quarks and different from 
and ±1 for exotic leptons. The systematic analysis of the 3-4-1 gauge theory carried out in fsl has shown that the 
restriction to fermion field representations with only ordinary electric charges allows for eight different anomaly free 
extensions. Four of them are three-family models and can be classified according to the values of the coefficients b 
and c which appear in the most general expression for the electric charge generator in SU{A)l <Xi U{1)x [see Eq. ([1]) 
below]. The allowed simultaneous values for these coefficients are b — c — 1 and b — 1, c — —2. Two of the four 
three-family models belong to the b = c = 1 class; the other two belong to the 6=1, c = —2 class. 

In 3-4-1 extensions without exotic electric charges the Z — Z' — Z" mixing can be constrained to occur between 
Z and Z' only, which leaves Z" = Z3 as a heavy mass eigenstate The diagonalization of the Z — Z' mass 

matrix produces a light mass eigenstate Z\ which can be identified as the neutral gauge boson of the SM, and a 
heavy Z^- After the breakdown of the 3-4-1 symmetry down to SU(3)c ® V(\)q, and since we have one family of 
quarks transforming differently from the other two under the gauge group, one important difference between the 
aforementioned two classes of 3-4-1 extensions appears: even thought in both classes the Z\ current remains flavor 
diagonal, in the b = c — \ class the new Zi gauge boson couples nonuniversally to ordinary quarks thus transmitting 
tree-level FCNC at low energies, while the couplings to Z3 are universal and couples only to exotic fermions. In the 
& = 1, c = —2 class, instead, it is the new gauge boson the responsible for this effect because couples nonuniversally 
to ordinary quarks, while the Z^ current remains flavor diagonal and universal. 
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The study of the bounds on the Z2 and masses coming from electroweak precision data and from FCNC processes 
has shown that extensions for which 6=1, c = — 2 are preferred in the sense that they give lower bounds on these 
masses smaller than the bounds predicted by extensions in the b = c = 1 class [9]. In fact, in the latter class a fit 
to Z-pole observables and atomic parity violation data produces m^^ ^ 2 TeV, a bound that is further increased to 
"^^2 TeV when the constraints coming from neutral meson mixing in the study of the tree level FCNC effects 

are taken into account. In the 6 = 1, c = —2 class, instead, the same fit gives the lower bound toz^ > 0.8 TeV and the 
analysis of the constraints arising from neutral meson mixing provides mz^ ^ 6.5 TeV. This means that extensions 
in the 6 = 1, c = — 2 class have a better chance to be tested at the LHC facility or further at the ILC. Here however 
we must point out that, as it has been noted repeatedly in the literature, the bounds from neutral meson mixing are 
always clouded by the lack of knowledge of the entries in the quark mixing matrices and involved in the new 
physics contribution, which forces to adopt well motivated but ad hoc textures for these matrices. 

In this work, in the context of a 3-4-1 extension belonging to the 6=1, c = — 2 class (the so-called "Model 
F" in Ref. 0), we re-examine the issue of tree level FCNC transmitted by the new Z3 gauge bosons but, due to 
the uncertainties associated with our ignorance of the quark mixing matrices, we will not search for bounds on the 
mass, but rather we set this mass at fixed values in the range 1 — 5 TeV (to be justified below) and search for 
information about the size of angles and phases in the Vj^ mixing matrix. To this purpose we will use several well 
measured AF = 2 {F — S,B) observables in the down sector, namely AAIk, AAI^i/s, £k, and sin$rf, associated to 
the Kl — Ks and the B'^^^ — 5°/^ mass differences, the Kaon CP-violation parameter and the — B^ mixing phase, 
respectively. We will also study the implications of the obtained bounds for the modifications in the amplitudes of 
the clean rare decays ir^vv, Kl — > 7r'^;//>, -B^/s ^ M^M^ ^^id Kl tt'^Z+Z^ {I = e,/i). In this context our 

mail goal will be to obtain upper and lower bounds for the corresponding branching ratios (BR) and to study their 
compatibility both with the experimental data and the SM predictions. 

Another well known FCNC rare process that deserves attention due to its sensitivity to new physics, is the radiative 
decay 6 — >■ 57. In general, this decay receives contributions both from the new charged and neutral gauge bosons 
and from the new scalar fields and can be used to put limits on these sectors of the 3-4-1 extension. As can be seen 
from the 3-4-1 scalar structure in Eq. below, and as it is done in the 3-3-1 model, the scalar contributions can 
be accounted for by an effective two-Higgs-doublet model. This, however, demands the identification of the physical 
scalar fields and their couplings and, therefore, the diagonalization of the full scalar sector. It is then clear that the 
study of this decay requires a particular and dedicated analysis. Moreover, as in the SM, 6 — > 57 is also a loop process 
in the 3-4-1 construction and since, as declared above, we are mainly interested in the analysis of the more stringent 
tree level FCNC processes that receive contributions from Z3 exchange, the detailed study of this decay in the context 
of the present SM extension will be postponed to a future work. 

This paper is organized as follows. In Sec. |IT1 we introduce the 3-4-1 extension to be considered and present its 
most important features for our purposes. In Sees. Illll and llVi we study the theoretical expression for the observables 
associated to neutral meson mixing and for the rare decays that will be analyzed, in such a way that the tree level 
FCNC contributions of the extra neutral gauge boson enter as corrections to the SM expressions. In Sec. |Vl we 
numerically evaluate the theoretical expressions obtained in the two previous sections, from which some information 
on structure of the mixing matrix can be obtained and the implications for the rare decays can be calculated. In 
the last Section we summarize our results. 



II. THE 3-4-1 EXTENSION 



Let us begin by introducing the most relevant features of the 3-4-1 extension we are interested in. Several details 
of their phenomenology have been already worked out in [l^l . 

The extension is based on the local SU{3)c <8) S'C/(4) l (8) U{l)x gauge symmetry which contains SU{3)c SU{2) l <8) 
U{1)y as a subgroup, and belongs to the 6=1, c = —2 class, where 6 and c are parameters appearing in the most 
general expression for the electric charge generator in SU{4)l <E) U{l)x 



Q = aTsL + -^bTsL + ^cTisL + X/4, (1) 

where TiL = XiL/2, with A^i, the Gell-Mann matrices for 5?7(4)l normalized as Tr(AiAj) = 25ij, I4 — diag(l, 1, 1, 1) 
is the diagonal 4x4 unit matrix, and a — 1 gives the usual isospin of the electroweak interaction. 

The anomaly- free fermion content of this extension has been discussed in Ref. [1] and is given in Table U where 
i — 1,2 and a = 1,2, 3 are generation indexes. The numbers inside brackets match the SU{3)c, SU{4:)]^ and U{l)x 
quantum numbers, respectively. Ui and C/3 are exotic quarks of electric charge 2/3, Di and are also exotic quarks 
with electric charge —1/3, while E~ and are exotic leptons with electric charges —1 and 0, respectively. 
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Table I: Anomaly-frcc fermion content. 



QiL ~ 

[3,4 


Ui 
T T 


d^L f/^i i^^L 

[3*,l,i] [3*,l,-i] [3*,l,-i] [3*, 1,^1 


QiL — 
[3, 


/ 1^3 \ 

dz 
Da 
\Ua) 


uIl d^L -Dsi UaL 
[3*,l,-i] [3*,l,i] [3*, 1,1] [3*,l,-i] 


LaL = 

[1,4 




[1,1,1] [1,1,1] 



The symmetry breaking occurs in three steps 



SU{i)c®SU{A)L® U{l)x 

SU{i)c®SU{'i)L®U{l)z 
A SU{i)c®SU{2)L®U{l)Y 

SU{i),®U{l)Q, (2) 

a task that is done by the following four Higgs fields with vacuum expectation values (VEV) aligned as 

(^D =((0;,0^0'i+,<^'i°)) = («,0,0,0)^ [l,4*,l/2], 
(0^) = ((02^,0^, 0!>",0r)) = (0,«',0,0) ^ [1,4*, -1/2] , 
Z-^^) ^§,^^0, = (0,0,^,0)^ [1,4*, -1/2], 



3 , 

) = ((00,0+, = (0,0, 0,n^ [1,4*, 1/2]. (3) 



This scalar structure consistently gives masses for all the gauge bosons and it is also enough to produce the observed 
fermion mass spectrum for the charged fermion sector (quarks and leptons) provided the hierarchy V, V' » v, v' ~ 
174 GeV is satisfied, where V and V set the mass scales for exotic fields [10| . 

The gauge couplings 54 and gx^ associated with the groups 5'J7(4)i and J7(l)x, respectively, are defined through 
the covariant derivative for 4-plets as: iD^^ = id^ — g^XLaA'^/'^ — gxXB^^. When the 3-4-1 symmetry is broken to 
the SM, we get the gauge matching conditions 

111 

34 = 3, and = — + ^, 4) 

9'^ 9x 9^ 

where g and g' are the gauge coupling constants of the SU{2)l and U{1)y groups of the SM, respectively. 

Clearly, FCNC in this extension can arise from the mixing of ordinary and exotic fermions. These FCNC and 
violation of the unitarity of the CKM mixing matrix can be avoided by the introduction of the following Z2 discrete 
symmetry which assigns charges qz to the fields as [l3] 

gz(Qai,UQL,daL,-^aL,e+^,(?!)l,02) = 0, 

9z(f/aL,i^^L,£;+i,03,04) = 1, (5) 



where a = 1, 2, 3 is a family index as above. 
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The gauge invariance and the Z2 symmetry ahow for the fohowing Yukawa Lagrangian in the quark sector 

2 3 3 3 

i—\ Q — 1 a— 1 a— 1 

3 3 3 

+ E /^fa^oi] + QIlC\'\>1 E + E ^3orfoL 

a — 1 a— 1 a— 1 

3 3 

+</.4 E ^IoUIl + 03 E '^s'^^^l] + h.c, (6) 

a— 1 Q— 1 

where the h' s are Yukawa couphngs and C is the charge conjugation operator. From this Lagragian we get, for the 
up- and down-type quarks in the basis (ui, 1*2, U3, C/i, [/2, J/s) and {di, d2, d^, Di, D2, respectively, 6x6 block 
diagonal mass matrices of the form 



Muu = ( ^t' Ju ) and M,n = ( ^^^^ Jo ] , (7) 



where, for V ^ V' and v ^ v' , the entries in the submatrices are 

^ hl^v and M^^ h%V- 



For the charged leptons we get the Yukawa terms 

3 3 



Mt.^ht.v and M^.c^h^V. (8) 



= E E LlLC[4>2hlpe+^ + cl^sh^pE;^] + h.c. (9) 

a=l 13=1 



From this equation we find a block diagonal mass matrix in the basis (ei, 62, 63, Ei, E2, E^) given by 

where now the entries in the submatrices are 

M^^^ = hl^v' and Mf^^/zf^F. (11) 

The mass matrices in ([7|) and (fTO|) exhibit a simple mass splitting between ordinary and exotic charged fermions, and 
show that all the charged fermions in this extension acquire masses at the three level. Clearly, by a judicious tuning 
of the Yukawa couplings a consistent mass spectrum in the ordinary charged sector can be obtained. In the exotic 
charged sector all fermions acquire masses at the scale V ^ V ^ vew = 174 GeV. We also remark that the tensor 
product form of the mass matrices M^u and Mdo in Q implies that they are diagonalized by unitary matrices which 
are themselves tensor products of unitary matrices. So, the discrete Z2 symmetry also avoids violation of unitarity 
of the CKM mixing matrix. The neutral leptons i'^^ and N^^ {a — 1,2,3) remain massless as far as we use only 
the original fermion fields shown in Table HI However, their masses and mixing can be implemented by introducing 
one or more Weyl singlet states „ ^ [1, 1, 0], n = 1, 2, without violating our assumptions, neither the anomaly 
constraint relations, because singlets with no X-charges are as good as not being present as far as anomaly cancellation 
is concerned. 

A look at the Yukawa Lagrangians in Eqs. ([6]) and ([9]) shows that each flavor couples to more than one Higgs 
4-plet and, consequently, scalar mediated FCNC also arise. However, since all the processes we will consider in this 
work involve external light quarks and leptons, these contributions are suppressed by the small Yukawa couplings 
associated to these fields. Notice that, for these processes, this will also be the case even in the absence of the 
discrete Z2 symmetry introduced in Eq. ([5]) because the Yukawa couplings will now be suppressed by small mixing 
angles. In the charged boson sector there are SM-like gauge bosons that do not mix with the other charged 
bosons and acquire a squared mass = {g\/2){y'^ + v'^) so that, with M\y = 80.399 ± 0.023 GeV [ll|, we have 

_|_ y/2 ^ — GeV. The ten remaining physical charged gauge bosons, namely: K^, V^, Y^, X^{X'^), 
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and K'^{K'°), acquire masses at the large scale V' ^ V and, at tree level, couple always to at least one exotic fermion 
[To| . This means that, for low energy processes involving ordinary fermions, the contribution of the new charged 
gauge bosons will be present only at loop level. Thus, we expect that the tree level FCNC effects transmitted by the 
exotic neutral gauge bosons will dominate. 

Our main interest, therefore, is in the neutral gauge boson sector which consists of four physical fields: the massless 
photon Ap, and three massive gauge bosons which come from the diagonalization of the mixing Zp, — Z'^ — Z'^. In 
terms of the electroweak basis, they are given by 

Af^ = SwA^ 



z'^ = -L(l-^^)V2(^^._2^)-^^i3^ 

Z"^ = 2AyV6 + A^,^/V3. (12) 

Since we are interested in the low energy phenomenology, we can choose V V' . If we also take ~ u', the 
current Z"^ = Z^ decouples from the other two and acquires a squared mass 

Mh = iglf^W^+v^) UM- This fact 
produces an enormous simplification in the study of the low energy deviations of the Z couplings to the SM families 
which now come from the diagonalization of the mass matrix 



+Cw 


Tw 

.7f 






+ (1 






CwA^^ 












—Sw 


Tw 
.V3 






+ (1 







2 / Sv^Sw 



where 5 — gx /<?4, and Sw = 5/\j26'^ + 1 and Cw are the sine and cosine of the electroweak mixing angle, respectively. 
The corresponding mass eigenstates are: = Z^^ cos 6 + Z'^ sin^? and Zl^ = —Z^ sin^? + Z'^ cos 9, where the mixing 
angle 9 is given by 



'^^^''^ - .^(l + 5g,)^ + y^C^-2.^ ' (1^) 



with C2W — Cw ~ '^w- Since V » v this mixing angle is expected to be very small. 

A. Neutral Currents 

The Lagrangian for the neutral currents J^(EM), J^(Z), J^(Z'), and J^(Z"), in the basis in which all the fields 
are gauge eigenstates, is given by 

-^'^C ^ eAMj^(EM) + -^Z^J^(Z)+.gxZ'^J4Z') 

Cw 

' Z"J^{Z"), (15) 



2\/2 



where e = gSw — 9x Cw V 1 ~ > 0. Calling qf the electric charge of the fermion / in units of e, the currents are 
given by 
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2 _ 2 

J^(EM) = -[u3J^,U3 + U3-f^,U3+^{u^-f^.U^ + U^-f^U^)] 
1 _ _ 2 _ 

i=l 

3 

a=l 

= Y.vHf. (16) 

/ 



J^.{Z) = J^,,L{Z)-SlrJ^{^M) 

1 - ^ - _ 



i=l 

3 



Q = l 



2 

+ U3Ll)j.U3L - '^{T'^diLliidiL - T'^UiL^^UiL 

1=1 



a=l 

-Tle-^l,e-^ - E-^1,E-^ + N^^lA)] " J^(EM), (18) 

and 

2 
i=l 

3 

+ ^(^>aL7/.i^aL + e^^l ^.e^^ - KlI^^Kl - ^aL7/.^°L), (19) 
a=l 

Eq. P7|) shows that we can identify as the neutral gauge boson of the SM because J^{Z) is just the generahzation 
of the SM neutral current. Moreover, from it is straightforward to see that the neutral gauge boson Z'^ does not 
transmit FCNC at low energy since it couples diagonally and universally to ordinary fermions. 

Notice, on the other hand, that J^{Z") is a pure left-handed current and that, notwithstanding the neutral gauge 
boson Z'^ = Z!^ does not mix neither with nor with Z'^ (in the case V' V and v' ~ v), it still couples non- 
universally to ordinary fermions. As a matter of fact, even though the Z^' couplings are diagonal, its couplings to 
the third family of quarks are different from the ones to the first two families. Thus, at low energy, we will have 
tree-level FCNC transmitted by Z'^ which are induced by fermion mixing. This means that in the corresponding low 
energy effective Lagrangian the chiral Z^' couplings will in general explicitly depend on the entries of the unitary 

matrices V^^^ that diagonalize the quark Yukawa matrices. So, in the analysis of these FCNC effects a convenient 

parametrization of the matrices must be chosen. 

At this point we must notice that if the mass splitting between Z2 and Zt, is so large as the constraints from neutral 
meson mixing suggest, we must be aware about the fact that the loop contributions coming from the exchange of 
Z2 can compete with the tree-level effects transmitted by Z3. But, as already discussed, these constraints are rather 
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unreliable, and since the mass matrix in (jl3l) can be exactly diagonalized we can estimate this mass splitting in the 
case we are considering, that is y ~ y >> w' ~ w which, in turn, implies a small mixing. So, neglecting the mixing 
we have: w gWV^ /T^ = 5l^V(l - T^) which, compared with = igl/2){V^ + v^) w {gl/2)V^, shows 

that in our approximation these masses are of the same order V and, consequently, tree-level FCNC mediated by Z3 
will dominate. 

B. The Effective Lagrangian 

We will use the formalism developed in Ref. [l^] where general expressions for calculating FCNC effects in models 
predicting the existence of one extra neutral gauge boson are presented. We start by generalizing this formalism to 
the case of — 1 extra neutral gauge bosons and then we restrict ourselves to the case — 1 = 2. 

The Lagrangian for neutral currents in Eq. (jlSp can be rewritten and generalized as 

£nc = -eA^J^(EM) - 31^1% J(i)^ ~ J2 (20) 

a=2 

where = Z denotes the neutral gauge boson of the SM and Z° are the new heavy Z bosons (which in general mix 
withes'). 

Following the notation of Ref. pJi] , the currents can be written as 



= EEv^^7. 4T'pl + 4C'ph (21) 



where the sum extends over all quarks and leptons and Pr^l = (1 ± 75)/2- ^^i l-. ~ '^R.L{i)5ij denotes the SM 

chiral couplings and e^^™! {m 7^ 1) denotes the chiral couplings of the heavy gauge bosons. 

After the spontaneous symmetry breaking the physical massive bosons Za are linear combinations of the gauge 
eigenstates 

AT 

Za — E UapZ'^, (22) 
/3=1 

where U is an orthogonal N x N matrix. 

The chiral Z^ couplings to fcrmions in the fermion mass eigenstate basis read 

i^t^^-<«45V,V, (23) 

where the CKM mixing matrix is given by the combination 

l^CKM = Vi:Vj^\ (24) 
The four-fermion effective Lagrangian, for the general case of A^ neutral gauge bosons, reads [13] 

--^--^E^MEt^-^-^^M ' (25) 

a=l \f}=l 31 

where pa = {m\cos^9w) and nia is the mass of Za- 

Using Eqs. (|211) and (l23l) into (|25l) a general expression for the effective Lagrangian, written in a way useful for 
further calculations, is obtained as 

- £eff = ^ E E E E {^a,Px^,) {XmrPvXn) , (26) 
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where k,l = 1,2,3; X and Y run over the chiralities L,R; ip and x represent classes of fermions with the same SM 
quantum numbers, i.e. u, d, e~ , v, while i, j, m, n are family indexes, ly^^^"^"- is given by 



W 



kl,ijmn 9k9l 



XY 



9k9l f ^ rr rr \ pi>{k) pX(0 

^1 \r=l / 



(27) 



where, in order to have a compact expression for the Lagrangian in Eq. (j26p . the summation is written so that takes 
elements of the matrix U. 

Another way of writing Eq. (|26l). using the notation of Rcf. 12], is 



-Ceff = 



4Gf 



mn 



(28) 



with the local operators given by 



mn 



mn 



Qi3 



{■iPiY'PriIJj) {Xmlt^PRXn) , 
{'tpij'^PLlpj) {Xmlp^PFtXn) , 
((ipil^PRtpj) {Xmlp^PLXn) , 



(29) 



and the coefficients are 



mn 



kl,ijmn 
LL ' 



kl,ijmn 
LR ' 



kl 



^mn — E ^RR 



tjmn 



kl 



mn 



kl^ijmn 
RL 



(30) 



The former expressions are perfectly general for any number of new Z bosons. In the case of only one extra Z 
boson, from the these expressions and using an orthogonal 2x2 transformation matrix U parametrized by a mixing 
angle 9, it is straightforward to obtain the formulae for C^^ , C^-J , D^^^ and D^^^ in Ref. fl^j . 



1. The N ^3 case 



Let us now restrict ourselves to the case TV — 1 = 2 corresponding to the 3-4-1 extension so that in Eq. (|20p we have 
Z2 = Z' , Z^ = Z" , and a comparison with Eq. ([T5|) gives: gi = g^/Cw, 32 = gx and 33 — 174/(2-^2). These gauge 
couplings are expressed in terms of the SU{2)l ® U(1)y couplings of the SM by the gauge matching conditions in 
Eq. (|4]) and, as it is done in studies of low energy FCNC effects associated to extra neutral gauge bosons when only 
SM particles are present in the effective theory, we neglect their renormalization group (RG) evolution from the high 
scale V ~V down to the Mw scale [see, for example, Refs. [13, El]]- 

Since in our case the Z — Z' ~ Z" mixing occurs between Z and Z' only, the matrix U takes the form 



/ cos e sin d 0\ 
U= -sine cose , (31) 
V 1/ 

Now, as already stated, the couplings of Z'^ to the third family of quarks are different from the ones to the first 
two families, thus allowing FCNC at tree level induced by fermion mixing. These FCNC have consequences on the 
predictions of AS* = 2 and AB = 2 observables in the down quark sector as the well measured AMk, AM^/s, sk 
and sin$d- In order to build the theoretical expressions both for these observables and for the related rare decays we 
will study let us extract, from Eq. p9)) . the couplings of Z'^ to ordinary down-type quarks. In the gauge eigenstate 
basis Q"^ — {did2d^)^ the interaction Lagrangian can be written as 
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E [^'7''(4S^^L + e%PR)Qf^ Z;, (32) 



where e^^"^-* = —-^3x3 + 2diag(0, 0,1) and e^^^^ = 0. Then, the chiral coupUngs, in the mass eigenstate basis 
Q'^ = {dsbY, are 

i^i^'' - -/3X3 + 2Vl diag(0, 0, 1)f/; Ef^ = 0. (33) 

Clearly, FCNC arise in this extension from the non-diagonal elements in the 3x3 matrix sf^^'^K That is, the 
coefficients which determine the contribution of £eff to AS* — 2 and AB — 2 processes are E'^*''^^ with i ^ j- 
Therefore, using ([25)1 and p3l) , the corresponding effective interaction Lagrangian is given by 

AGf 



^off = 



V2 \gi 

The nondiagonal elements of sf^'^'^ in Eq. (pS)) read 



^ A3 (Qf7.i?flf P.Q,^) (Ql7-<^i:PLQf.) . (34) 



- 2yL.3^Ll3, (35) 



whereas the ratio (33/51) can be written in terms of the Weinberg angle as 



91 



(36) 



As it is well known, in the SM the matrix entries in are not observables while the observable quantities only 
concern the entries in the Vckm matrix. However, from Eq. (|34|) it is clear that the new physics contributions will be 
proportional to the entries V^^^ in the down quark mixing matrix. We must, then, conveniently parametrize the 
matrix so that we can estimate the size of its entries. This, in turn, requires to determine the number of independent 
parameters in Vj^. An unitary 3x3 matrix has, in general, nine independent parameters: three real angles and six 
complex phases. In the CKM mixing matrix five phases can be absorbed in redefinitions of the quark fields, which 
leaves us with three real angles and one complex phase. In the case of the matrix, by analogous redefinitions of 
the three ordinary down-type quarks, we can absorb three complex phases so that we remain with three real angles 
and three complex phases. The matrix can consequently be conveniently parametrized as the product of three 
rotations each one of them containing a complex phase [ll|. We thus get 



-s°3e"*23 



5^4 

Performing the product we have 




= I c^3 _ ^236-^''- II 1 II _.sf,e'^'- cf^ I (37) 



„d d d d „ — ii5,, d 

--12'-13 *12'-13^ *13' 



yd \ d „d „d „d „<i „i(5f3-4,) „d „d „d „d „d piiS'is-S'^^-S^.^) „<i „d ^-iS^^ /'QQ^ 

" I -S;^2C236 - Ci2S23*136 ^'^ C12C23 - Sl2S23'Sl36 ^2301^6 23 | (^OOJ 

„d „d fiiiSfn+Sis) ^d „d „d AS'}„ „d „d „iS^„ „d „d „d „i{Sf^-S'^^) „d „d 

512*23^ ~ Ci2C23*13'S ~'^12^23^ ~ Sl2'^23^l3,^ '^23.'^13 



Another way Vj^ can be parametrized is |l7j 



Vi = PVK, (39) 



where P = diag(e'^i , 1, e*^^), K = diag(e'"% e*"^, e*"^) while, using the standard parametrization, the unitary matrix 
V can be written in terms of three mixing angles 612 , O23 and ^13 and one phase ip [ll| . Note that this parametrization 
includes six complex phases but again, by redefinitions of the ordinary down-type quarks, we can absorb three of them. 

It is easy to show that the theoretical expressions resulting from both parametrizations are equivalent via redefinition 
of phases [See Eqs. (|62II64I) below and Ref. 0]. 
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Table II: Values of the experimental and theoretical quantities used as input parameters. 



input 


\T 1 
VcllXlS 






L.lK)K)6(yl) X iU 




a 


T on'yQ coc o'7i5 / \y ia — 3 
( {6oZD6{b\p\}) X iU 




o^s (rnz) 






sm 6w {mz){MS) 


0.23116(13) 


LIU 


\Vud\ 


0.9728(30) 




\Vua\ 


0.2250(27) 




Vub| 


0.00427(38) 




\Vcd\ 


0.230(11) 




\Vcs\ 


0.98(10) 




\Vcb\ 


0.0415(7) 




n/ 

1 


A f\C\Q A ( ti\ 

U.UUo4(D j 






A flQQV/'OI \ 

U.Uoo / l^zl ) 








n 11 


Tow fnpvi 

/ / 1 VV 1 VIC V J 






TO 7 frieVl 


91 187fi('91~l 




TO [Hp VI 
I tic. ['-Jt; V J 






TO* rnpvi 


1 79 zl/i 9~\ 


[111 
LLLI 


7 


7sr 1 9^ 




J K V [vjC V J 






/s,\/Bi:[GeV] 


216(15) X 10"^ 




/s, /B^ [GeV] 


275(13) X 10"^ 


[20] 


»7i 


1.32(32) 


[22, 23] 


»72 


0.5765(65) 


[22, 23] 


»73 


0.47(4) 


[23, 24] 


»7b 


0.551(7) 


[22, 25] 


Kg 


0.92(1) 


[20] 



m^o [GeV] 497.614(24) x 10" 

ruBa [GeV] 5279.50(30) x 10"^ 

niB, [GeV] 5366.3(6) x 10"^ 

Amif [ps~^] 0.5292(9) x 10"^ 

Arns^ [ps-i] 0.507(5) 

Ams, [ps-^] 17.77(12) 

|£k| 2.228(11) X 10"^ 

sin$d 0.673(23) [11] 



III. AF = 2 OBSERVABLES 

We now proceed to build the theoretical expressions for the AF — 2 {F = S, B) neutral meson mixing observables 
AMk, AMd/s, ek and sin<l>d in such a way that the Z'^ contributions enter into the ejqDressions as corrections to the 
SM predictions, following a similar procedure as for the minimal 3-3-1 model in Ref. Q and for the economical 3-3-1 
model in Ref. Sj] . These expressions are functions of the matrix element 

where P stands for Bg or B^. In our case Mj^ receives both SM contributions arising from standard one loop 
diagrams and contributions coming from tree-level Z3 exchange, that is 

M[,^m[,'-^^^+m[,^'-'-'\ (41) 
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The expressions for the AF = 2 observables are 

ArriK = 2Re(Mj^), 
Anid = 2 

Am. = 2 M- 



^"12 
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\ ArriK 



Po 



arg(Mit), 



(42) 
(43) 

(44) 

(45) 
(46) 



where the term Pq is due to long distance contributions to Kaon mixing. The experimental values for these observables 
are consigned in Table HIl It is customary to set (pe = vr/4 and to neglect Pq flo, To']. For this reason a multiplicative 
correction factor for ek, that accounts for ^ 7r/4 and Pq 7^ 0, is introduced (2^ i21(]. Then, the expression for ek 
becomes 



Ke exp 



f.n 

V 4 



7r\ Im(M- 



/ ^/2A 



rriK 



The well known SM contributions to M/^ are given by [3 [l^ 



(47) 



i<:(SM) 



r 

TTT^fKBKmKm^ A* 7yi5o(xc) + A^* 772S'o(xt) + 2A*At 773S'o(a;c, Xt) 



and 



^B,(SM) 



^ fB,BB,VBmB^m'^So{xt){Vt*gVtb)'^, 



127r^ 



(48) 



(49) 



where g stands for d or s. /p is the P-meson decay constant, Bp is the corresponding renormalization scale and scheme 
invariant bag parameter and the basic electroweak loop contributions, without QCD corrections, are expressed through 
the functions So{xi) {xi — m^/m^), namely 



S'o(xc) = Xc, 



(50) 



So{xt) 



Axt 



llxi 



Zx^ In Xt 
2{l-xtf' 



(51) 



So{Xc, Xt) = Xc 



In 



Xt 



ixt 



Sec? In Xt 



4(1 -Xt) ^{l-xtf 



(52) 



Renormalization group short-distance QCD effects are described through the renormalization scheme independent 
factors ?7i, 772, % and -qB [22| - [25j . and the CKM factors are given by A^ = V*.V^^. 

Now, from the effective Lagrangian in Eq. ([M)) the expressions for the 3-4-1 contributions can be obtained. Since 
only left-left (LL) operators appear in the effective Lagrangian, the QCD RG evolution of the Wilson coefficients C(/i) 
after including the Z3 contributions is exactly the same as in the SM. Then, the relevant elements of the evolution 
matrix U{n, /X3) from the high scale /X3 — mz^ down to the scale /Lt, in the notation of Ref. [26|, can be calculated from 



[?7(m,M3)]vll = [r|[^l,^lt)\^rh\.Vl[^J't,^J■3)\ 



VLL, 



(53) 



where /ij — 0{mt). 

In general, for /i < /i', the correction factor [77(/i, /i')]vLL in an effective theory with / active quark flavors is given 

by 



„(/) 



[r?(M,/^')]vLL = [r?^"HM)]vLL + 



47r 



VLL, 



(54) 
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where [??^°Hm)]vll and [?y''^HM)]vLL are the leading order (LO) and the next-to-leading order (NLO) factors, respec- 
tively, and in the right-hand side we have suppressed the high scale /i'. 

For K and B meson mixing the factor [r]{iJ,, Ht)]vLL in Eq. (f53|) has been calculated in Ref. [IHl using = mt{mt) = 
166 GeV. The results are: [?7^(AtL, /it)]vLL — 0.788, for — K'^ mixing with fi^ — 2 GeV the lattice scale where the 
matching with the lattice results of the associated hadronic matrix elements is done, and [f7'^(/i6, /it)]vLL = 0.842, for 
^d/s ^ ^d/s mixing with fib = 4.4 GeV the mass scale of the bottom quark. In order to calculate [77(^x4, /is)] vll we 

first evaluate 0^^113) to LO in the MS scheme in an effective theory with f — 6 quark flavors and we impose the 
continuity relation af \fit) — ai^''(/it). To this purpose we use 

1 - 1 >f^o\^,t^^ 



In^, (55) 



where = 11 - (2/3)/ and for a'fUmz) we use the central value ai^\niz) = 0.1184. In this way [?7(/it, /i3)]vLL 
can be obtained from Eq. ([Ml) with [2g 



VLL 



^^'HmOJ^^^ = . (56) 

where r]G = aT\fi3)/ah\fit)- The numerical value of [77(/it, a*3)]vll depends of the value of the high energy scale 
M3 = TUz^- Since in our approximation the Z2 and Z3 mass scales are of the same order (see Sec. |II]), we will select 
values for mz^ in the typical range 1 — 5 TeV. 
With this in mind and defining 

C^VLL = [?7^(ML, )] VLL [?y(Mt, Ma)] VLL, 

C^VLL = [?y^(Mh,Mt)]vLL[r?(Mt,M3)]vLL, (57) 

the 3-4-1 contribution to Af/^ is given by 

,,fP(3-4-l) SV^Gfj^p f 93\^ .2d \2 rro\ 

*^i2 = o t^VLLPsI — ) mpfpBpXp, (58) 

where 



3 \gi 



Xk = 4 4 e^* , (59) 
As. = 43 4,4,e^'^', (60) 
As. = 44(c^3)'e-^^^ (61) 



and the phases are defined as 



= -5f„ (62) 

0" - -(5^3 -^23), (63) 
03 = 0" -0'. (64) 

As can be seen from Eq. ([55)1 . the 3-4-1 contributions M^^^^ are given in terms of five unknown independent 
parameters, namely, the mass of the gauge boson, the angles ^^3, 6*23 and the two complex phases (j)' and (f>" , where 
the last four parameters come from the Vf mixing matrix. 

The AF — 2 observables discussed in this Section are, in principle, sufficient to set bounds on the relevant parameters 
of Vl in Eq. 
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Figure 1: Allowed S23 — 4>3 region for sfs — 0, and mz^ — 1 TeV (Left), mz^ — 5 TeV (right). 



Table III: Upper limit for 6\z and 62?, for different values of mz^- 



mz3 


^13, max 


^23, max 


1 TeV 


8.73 X 10"" 


3.14 X 10~^ 


3 TeV 


2.97 X 10"'' 


9.42 X 10"^ 


5 TeV 


5.06 X 10"^ 


1.75 X 10-' 



IV. THEORETICAL EXPRESSIONS FOR RARE DECAY AMPLITUDES 

In the next Section, using the observables related to meson mixing, the bounds on the parameter space associated 
to the relevant entries of the down quark mixing matrix will be established. Our next task will be to study the 
implications of these bounds on several rare decay amplitudes following a similar approach as in Ref. Q . The decays 
we will study are the theoretically very clean rare decays — > tt^z/i?, Kl — tt^i^P, and the clean rare decays 
Bd/s ^ fJ-^fJ-^, Kl n"l+l- {I = e, fi). 

As it is known, the rare decays we are interested in are governed both by electroweak- and photon-penguins and 
by leptonic box diagram contributions. In the SM these contributions are described, in the LO, by the corresponding 
Inami-Lim functions Co{xt), Do{xt) and BQ{xt) which, in the expressions for the decay amplitudes, always appear in 
the gauge invariant combinations: Xoixt) = CQ{xt)-4Bo{xt), Yo{xt) = Coixt)-Bo{xt), Zo{xt) = Co{xt) + {l/4:)Do{xt) 
[27| . In the NLO these combinations will be written as X{xt), Y{xt) and Z{xt). Since the Z3 contribution to the 
effective Hamiltonians for the rare decays we are considering have the same operator structure than the ones in the 
SM, that is they have a {V — A){V — A) form, its effects on the various decays can be encoded by appropriate 
redefinitions of the X{xt), Y(xt) and Z[xt) functions in the form 

X'{xt) = X(^^'^)(a;,)+AX, 
Y'{xt) = y(^^'^)(xt) + AF 

Z'(xt) = Z(^*^^(xt)+AZ. (65) 

As we will see, in the 3-4-1 extension the redefined functions are in general complex. We will write them as they appear 
in Ref. fis'l so that the modifications lead to the correct results in the 3-4-1 extension including the corresponding 
imaginary part. For the decays — >■ n'^l^l^ , in which there is also a right-handed contribution, the Z3 effects will 
be absorbed into the matching conditions of the Wilson coefficients. 
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Figure 2: Allowed sf^ — (f)' region for mz^ 
0.002 (center) and 0.005 (right). 



1 TeV and for different values of (j)" and increasing values of sfg: 0.001 (left), 



For K — > ni^i' there exists a charged decay tt+i/i^, and a neutral one — > n'^i'iy [1^. For both of them the 

NLO effective Hamiltonian is given by 



ttSM Gf a 



y/2 27T sin^ ew 



J2 (KiKd^k + v,:vt,x{xt)) {-sd)v-A{i>iW, 



V-A , 



(66) 



where the functions Xnl and X{xt) represent the charm and top- loop contributions, respectively. Summing over the 
three neutrino flavors and collecting the charm contributions in PdX) — 0.41 ± 0.05 [23|, the BR for A'+ tt^vD 
can be expressed as 
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Figure 3: Upper and lower bounds on the BR for the decay — >■ Tv'^iyv (left), and for the decay Kl — >■ tv'^i'v (right). 
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Figure 4: BR(ii'+ -5> 11+ uu) /'QR(K+ -5> 7r+!/i5)g^ (left) and BR(ii'L -> 7r°!/!?)/BR(J<L 7r''i^!?)gj^ (right). Upper and lower 
bounds for different values of mz-, are shown. 



BR{K+ Ti+vv) = n+ 



Re ( -fPciX) 



+ Re ( -^Xixt) 



where |30| 

K+ = rK+ ■ 
For Kl tt^vv it is found 

with 



27r2 sin^ 



w 



= (5.26 ±0.06) • 10 



-11 



A 



0.225 



BR{Kl ^ 7r%9) = KL ■ (lm(^^X{xt)^^ , 



Kl = K+ 



rLc+ t{K^ 



(2.29 ±0.03) • 10" 



0.225 



(67) 



(68) 



(69) 



(70) 
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where A is a parameter of the Wolfenstein parametrization of the CKM matrix [3l| and is set equal to S12 of the 



standard parametrization 11|. 



For both decays the contribution from the 3-4-1 extension to the effective Hamihonian comes from a tree diagram 
transmitted by the Z" boson and has the form 

i/f;= ^f-Vp3l^/23l^/*13(^-d)y-A(^;^0v-A. (71) 

This contribution can be included as a modification AX of the function X{xt) which reads 

a \9iJ ytlVtd 
B. Bd/s -5> 

These decays are entirely determined by the top contribution so that the SM effective Hamiltonian in the NLO is 
given by 

= -^^iV;,V,,^,)Yixt)ibq)v^A{flf^)v-A . (73) 
From here, the expressions for the branching fractions are obtained as 



BRiB, ^ M+^-) = rs.^ms, f ^l^^V M^^^\V,IV,,Y ix,)\\ (74) 
The 3-4-1 contribution to -6^/^ M^M^ is found to be 



z" _ 2Gf ( 93 

'eff — 7^ 

which we now absorb into a modification of Y{xt) as 



2 

d T rd* 



H^S ^^[fj P3^L33^L13/23(^9)v-A(AiM)y-A, (75) 



"^T^Sm 9w f 93\ ^L33^L13/23 

AYb^^p, = — P3— TTTTT • (76) 



In the SM the CP- violating part of the effective Hamiltonian, after neglect QCD penguin operators, is given by 

H,^/'' = -^v*Vtdiy7vQ7V + ViaQia) , (77) 

where Q-jv = {sd)v -AG^'^e and Qta = (sd)v-AG''^^l^e. are the vector- and axial-vector operators which contribute 
and that originate from 7- and Z-penguin and box diagrams. The matching conditions of the Wilson coefficients y-jv 
and yjA are given by 

V7V = ^ - 4Zo(xO + Po) , (78) 



27r Vsin^ 0w J 
and 

a Ya{xt) . . 

VTA = -TT- . 2n ' ('^9) 

ZTT sm Uw 

In y-^v a small term Pe has been neglected and we use the normalization Pq — 2.89 ± 0.06 [32l. l33j. 
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For these decays the contributions from new physics lead to the effective Hamiltonian 

2 



%/2 V51 



(80) 



Following [j], instead of absorbing these new contributions into modifications of the Inami-Lim functions, they will 
be absorbed into the matching conditions of the Wilson coefficients in the form 



9A\ (^/23^/*13) 



(81) 
(82) 



The expressions for the BR, inclu ding the long-distance indirectly CP-violating terms and their interference with 
the short-distance contributions, are [SJ 



-"CPC 



where 



2 I , .2 \ 

lA) 



Qi, = (4.62 ± 0.24) (w^y-hw, 
C?t = (11.3±0.3)w7y, 
C^i^ = 14.5 ± 0.05 , 



c: 



CPC 



0, 



with wtajv = Im(At2/7yi,7y)/IniAf. 



10 



-12 



/'I nn _i_ n nc;^ .2 

CL - (2.63 ± 0.06)w7y , 



di,- (1.09±0.05)(u;7V 

int 



2.32w^^) , 



= 3.36 ± 0.20 , 
1.2 ±0.2, 



Cgpc = 5.2 ± 1.6, 



(83) 



(84) 
(85) 
(86) 
(87) 




BR(K")[xlO" 



Figure 5: BR(A'l -> ti^^vv) ~ BR(A'+ -> tt+i^i/) plane for mz^ = 1 TeV. 



V. CONSTRAINTS AND IMPLICATIONS FOR RARE DECAYS 



In this Section we numerically evaluate the theoretical expressions obtained in the two previous Sections. As 
noted above, the 3-4-1 contributions M^2 ^ ^ to the matrix elements in Eqs. ()42M6p are given in terms of the five 
unknown independent parameters mzs, the angles 6f^, 6*23 and two complex phases, </>' and cj)" . We perform two 
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Figure 6: Upper and lower bound on the BR for the decay Bg — > /i^/i (left), and for the decay Bd — >■ M^A* (right) as a 
function of mz^ ■ 
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related numerical analysis: in the first one, we use the well-measured observables AMk, AM^^/^, sk and sin$(j to 
constrain the angles and phases appearing in the Vj^ mixing matrix for selected values of the Z3 mass. In the second 
analysis we study the implications of these bounds on the rare decays previously mentioned. In this context we are 
mainly interested in obtaining upper and lower bounds for these decays as a function of the mass in order to 
confront them with the present experimental data and with the SIM predictions. The updated experimental data for 
the input parameters are also given in Table HIl 



A. Constraints from - K° and B° - B° mixing 

We start by focusing on the bounds on the parameter space of the 3-4-1 extension associated to the down-like quark 
mixing parameters. Two possible analyses can be carried out. The first one, followed in several works both in 3-3-1 
models [2] and in 3-4-1 extensions 0, 0, consist in to assume a texture for the mixing matrix in order to obtain 
bounds on the mass of the heavy gauge bosons. In the opposite approach, one can fix these masses and obtain some 
information on the size of the corresponding mixing matrix elements • Here we follow this second approach with 
fixed values of mz^ in the range 1 — 5 TeV. 

In previous studies on flavour physics observables in different models, the numerical analysis was simplified by 
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Table IV: Present experimental data and SM predictions for the rare decays considered in this work. 



Label 



Decay 



Experimental 

-1(1 



SM 



BR(_R'+) 

BR(B,) 
BR(Bd) 
BR(A'ie) 



BR{K+ TT+w) 
BR{Kl -> TT°iyu) 
BR(B, -> fi+fi-) 

BR(Bd ^ /i+M") 
BR(7^^L -> TT°e+e-] 



(1.7± 1.1) • 10 

< 6.7 ■ 10-** (90%CL) 

< 1 • 10"^ (90%CL) 

< 3 • 10^* (90%CL) 

< 28 • 10"" (90%CL) 



{l.Slton ±0.29) • 10 
{2A3t°oi°7 ± 0.06) ■ 10" 
(3.35 ± 0.32) • 10"'' 
(1.03 ±0.09) ■ 10^"^ 
(3.54tr8^).10-" 



< 38 • 10" 



(90%CL) 



(i.4r 



10" 



setting all input parameters to their central values and allowing instead AM^, e^, AM^, AM^, \sii\ and sin^^ to 
differ from their experimental values by ±50%, ±40%, ±40%, ±40%, ±20% and ±8%, respectively. This simplifying 
assumption was justified in order to determine the size of possible effects on observables that have not been well 
measured so far. In recent analysis, however, an improved error analysis has been done that enable us to draw more 
accurate conclusions in view of the recent significant improvements in the experimental constraints and in the lattice 
determinations of the non-perturbative parameter. Therefore, in what follows, we will take all input parameters to 
be flatly distributed within their la ranges indicated in Table HIl At the same time, we require the observables \sk\, 
AMd, AMs and sin^^, resulting from the SM and the 3-4-1 contributions, to lie within their experimental la ranges. 
In the case of AMk, where the theoretical uncertainty is large due to unknown long-distance contributions, we allow 
the generated value to lie within ±30% of its experimental central value. 

In dealing with the values of the CKM matrix parameters quoted in Table |lll instead of using the full CKM angle 
fits, we take only the values related to direct experimental measures with the largest uncertainties taken into account. 
In this way we avoid possible modifications from the 3-4-1 contributions to the one- loop SM processes included in the 
global fits and leave the largest room available for possible new physics effects, respecting the well measured values. 
For the CP-violating parameter 7 we take its direct determination from the model-independent UTfit analysis of 
B decays. 

We then proceed to scan the parameter space fixing mz^ and generating a large number of points that we call an 
"effective" parameter space. For each fixed value of the Z3 mass, this space consists of more than 1 x 10^ points that 
fulfil the requirements listed above, where all angles are varied in the interval [0,7r/2], all phases between and 27r, 
and all input parameters are varied in their la ranges. 

With these data we can plot contours setting bounds on some of the parameters. For example, let us take the 
mixing angle = 0. In this case, Eq. (j6ip determines the allowed region in the S23 ~ 'Ps plane, which is shown in 
Fig. [1] for 771^3 — 1 TeV and mz^ — 5 TeV. These plots allows us to set upper bounds on jsjal which depend on the 
value of mzs- It is seen that the upper bound for S23 increase for increasing values of mz^- Fig. [T]also shows that, 
at the la confidence level, there is a large and typical region which is excluded and that prefers values for 03 around 
48° and 136° independently of the Z3 mass. 

Now we explore, for several values of mz^, the allowed regions in the plane sf^ — (j)' using non-zero values of S23 
and selected values of cj)" . This is shown in Fig. [5] for mz^ = 1. We have taken values for S23 consistent with the ones 
in the allowed S23 — (t>^ plane and we have fixed (p" at the values 0°, 45°, 90° and 135°. This figure shows that for 
increasing values of S23 and depending of the (j)" value, the allowed region reduces considerably involving both the 
selection of very specific values of 0' and a decrease in the upper bound for 5^3. 

With the data from the effective parameter space we can estimate the order of magnitude of the upper limits for Ofr^ 
and 6'23 for several values of the Z3 mass. The results are collected in Table [TTIl which shows that these upper bounds 
increase for increasing values of mz^ ■ This allow us to elucidate the size of the allowed region in the parameter space. 



B. Implications for Rare Decays 

In order to study the implications of the obtained bounds for the modification in the rare decay amplitudes given 
in Section llVi we move over all the effective parameter space for each selected value of mz^ , and we calculate the 
amplitudes for the different decays. At difference with the analysis in Ref. [3], our procedure allows us to set not only 
upper bounds but also lower bounds on the corresponding BR as a function of the Z3 mass. The present experimental 
data for these decays [ll| and the SM predictions [2^, [SJ, HI] are given in Table IIVI Our goal will be to determine 
these upper and lower bounds and establish their compatibility with the data in this Table. 

For the decays tt^vv and Kl ^ tt^vv^ using Eqs. (|57)) . ([M)) . ([7^ and the effective parameter space, the 
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Figure 8: Upper and lower bound on the BR for the decay Kl — > Tr^e^e (left), and for the decay Kl — )■ n'^fi^fi (right). 
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Figure 9: Upper and lower bounds for the ratios ^H{Kl — > ■K'^e^e )/'SR,{Kl yr^e^e )gj^ (left) and 
BR(A'i, — >■ tt'' )/BR(if I, — )■ 7r"/^"'"/i~ )gj^ (right), for different values of mzg. 



bounds are given in Fig. |31 For the decay K'^ — ^ t:^vv we see from Table HVl that the experimental central value is 
greater than the SIvI prediction. Fig. |3] shows that for low values of mz^ the upper limit reaches the experimental 
central value, but for larger rn^a both the upper and lower limits go closer to the SM prediction. For the decay 
Kl 11^ vv we can see that big enhancements are expected for low m^,. The bigger enhancements are expected for 
values less that ~ 2.5 TeV. For the decay — >■ n'^vv the greatest contributions come from the CP-conserving 
case and for Kl — > n'^vv come from the CP- violating case in accordance with the setting in 4], but with the difference 
that we have used all the effective parameter space which includes the CP-conserving, the CP-violating and the mixed 
cases. 

To appreciate more clearly the departures from the SM, the ratios BR(_ft'+ — > tt^vv) /'SK{K^ ^^'^'^)sm ^^"^ 
BR(i4'i — >■ 7r°i/P)/BR(i^i — 'K^vv)^^ vs. mz.^ are shown in Fig.|4l For the decay A"+ tt^vv we find enhancements 
around 2.8 times the SIM prediction for low values of mz^, but small departures from the SM for large mass. For 
Kl —7- TT^vv, enhancements of around 7.6 times the SM value are found for low values of m^g. 

Comparing the upper and lower bounds on both decays we see that the possibility of enhancement is always larger 
than the possibility of suppression, except for large m^, where both of them are comparable. 

The plane BR(i^i — )• tt^vv) — BR(i4'+ — ^ tt^vv) is shown in Fig. [5j The upper and lower branches correspond 
mainly to the CP-violating case and the horizontal branch to the CP-conserving case. It is implicit in this figure that 
the the mixed case is very restricted making the expressions in Eqs. (j72p . (j76|) . (1811) and (1821) mainly real or imaginary. 
Notice that the pattern seen is similar to the one obtained in the Littlest Higgs model with T-parity in Ref. [3gl and 
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Figure 10: BR{Kl -5> n° ^^+ ^i' )-BR{K l -5> vr^e+e") plane for = 1 TeV. 



in the minimal 3-3-1 model in Ref. f3|. 

For the decays Bg /^^M^ ^i^d Bd using Eqs. (|74|) and (|76|) . the limits are shown in Fig. [6l From the 

data in Table HVl we see that, for all the selected range of values of m^,, the upper bounds are smaller than the 
experimental upper limit. For the decay Bd — ?> fJ-'^ fi~ the upper and lower bounds are lower than the SM predictions. 

BR{Bs — > /i+^^)/BR(i3s — > A*^M~)sM and BR(i3rf — > fi^^^)/BR{Bd — > M^A* )sm ™-Z3 are shown in Fig. [71 For 
Bs — > fJ-^ fJ-^ we see that both enhancements and suppressions are possible but both of them decrease with increasing 
values of mz^ and are very small with respect to the SM prediction for mza = 5 TeV. For Bd — fJ-^ fJ-^ suppression of 
around 0.87 times the SM value is possible for mz^ ~ 1 TeV, but for higher mzs the suppression can each 0.92 times 
the SM prediction. 

Going back to Fig.|6]it is interesting to note that, due to the small difference between the upper and lower bounds 
for BR(i?s -> /i^/i^) and JiK{Bd — ?> ji^ jj."), they are suppressed in this 3-4-1 extension even with the inclusion of new 
CP-violating phases, as compared with the case for the decays in Fig. [3l A similar behavior is reported in [4] for the 
minimal 3-3-1 model. It is important to remark that, as in the 3-3-1 model, a strong enhancement of this BR would 
rule out the extension we are studying. 

We conclude analyzing the results for the decays Kl — > 7r°e+e~ and Kl -> 7r°/i+/i^. Using Eqs. ((8T|) . ((82)) and 
, the resulting limits are shown in Fig. [S] At this point it is important to clarify that the upper bound is built 
using the -I- sign in Eqs. ([55)) and the lower bound is obtained using the — sign. From the data in Table [TV) we can 
see that these limits are lower than the SM predictions, but consistent with the upper experimental bounds. Notice 
also that the upper and lower limits are both almost independent of the Z3 mass. 

BR(Xi 7r°e+e")/BR(ifi -)■ 7r°e+e~)gj^ and BR{Kl 7r°/x+^")/BR(i4ri n"fi+^i^)g-^ vs. mzs are shown in 
Fig.[Slfrom which we can see that these decays are suppressed compared to the SM prediction. Moreover, a comparison 
with Fig. [3 shows that, in the 3-4-1 extension, these BR are more suppressed than the ones for the Bg — > and 
Bd — > M^M" decays. 

In the spirit of [Hi, we plot the plane BR{Kl tt" ^1+ pi-)-BR{Kl Ti^e+er) for mz, = 1 TeV in Fig. [TUj This 
shows the correlation between these decays. It can be seen that they are very restricted. Comparing Fig. [TUl with the 
corresponding figures in Refs. ,32, 34], it is noted that they are very close to the plots for the SM in these works. This 
fact can also be observed from a comparison between the two plots in Fig. [8] 

In Fig.[IIl the planes BR(i^L ^ i:°e+e-)-BR{KL n^fp) and BR{Kl n° n+ fi-)-BR{KL n°iyP) are shown. 
We find that, as it happens in many other models, the decays — >■ Tr^/x^/x" and Kl — >■ 7r°e~''e~ are subject to weaker 
modifications than the decay — >■ tt^vP. We remark that this is contrary to the result found in Q for the minimal 
3-3-1 model. This result, then, can allow to test the 3-4-1 extension studied here when more accurate measurements 
of both decays become available. 
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Figure 11: BR(A'z, -> 7r°e+e-)-BR(A'z, ir^uiy) plane (left) and BR(i4:i, tv° ii+ )-BR{K l 7r°!/i5) plane (right) for 
mz^ = 1 TeV. 



VI. SUMMARY 



In this work, in the framework of a 3-4-1 extension of the SM characterized by the values b — 1, c = — 2in the 
most general expression for the electric charge generator and with chiral anomalies cancelling between the families, 
we have analyzed the effects of FCNC transmitted by one of the two extra neutral gauge bosons predicted by this 
construction in the case in which this new gauge boson does not mix with the other two neutral gauge bosons. In 
this context we have studied in detail the constraints that well measured AF = 2 observables {F = S, B) impose on 
the parameter space associated to the Vj^ matrix describing mixing in the down quark sector. This parameter space 
consists of five variables, namely the angles 6*^3, ^23 ^^.d the CP-violating phases </>', </>", related to the mixing 
matrix, and the mass, and we have used the well measured quantities IS.Mk, AM^j/^, £k and sin$d- For fixed 
values of some of the relevant parameters, allowed regions for the remaining parameters were plotted in Figs. [T] 
and [2] when mz^ is allowed to vary in the typical range 1 — 5 TeV. A calculation of the upper limits for ^13 and ^23, 
for different values of mz^ , is presented in Table IIIII which shows that these bounds increase when increases. 
This allowed us to appreciate the behavior of the parameter space for selected values of the Z^ mass. 

We have then taken these results to obtain upper and lower bounds for the BR of several clean rare decay processes. 



namely if"*" 



Kr 



Kr 



{I = e, fi) and B, 



d/s 



We find sizeable enhancements 



and/or suppressions with respect to the SM prediction mainly for mz^ values in the range 1 — 2 TeV. In this range 
large enhancements, as compared with the suppressions, are predicted for the BR of the decays — > tt^Di' and 



Kl -5> 7r°z/P. For BR{Bs 



) we find enhancements slightly larger than suppressions, while for BR{Bd — > /i ) 



we find suppression for all the considered range of values of mzg. Finally, for BR{Kl — >■ Tr^l^l~) we find upper and 
lower bounds smaller than the SM prediction but consistent with the upper experimental limit. These departures 
from the SM, presented in Figs. SI [7]and|ni can be considered as signals of the 3-4-1 extension under consideration 
when looking at the data. From these plots is clear that the decays -k^vv and Kl n^vv show the greatest 

differences with the SM predictions, but the decays Kl — >■ 7r°l+Z~ {I = e, n) and fi^/s M"*"/^" are very restricted in 
the 3-4-1 extension. It is important to note that, as in the minimal 3-3-1 model, a strong enhancement of the BR for 
the decays B^/s — > M^A*^ would rule out the 3-4-1 extension studied here. 

The plane BR(/Ci — >■ n'^i/p—) — BR(_fir+ — )• tt'^Di') in Fig. [5] allowed us to compare our results with the ones in the 
Littlest Higgs model with T-parity [36|] and with the ones in the minimal 3-3-1 model [3], and it also helped us to 
better understand the structure of the expressions in Eqs. (f72|) . (|76|) . (|8T|) and ([82|) . It is important to remark that, 
contrary to what it is done in [4] , our results are obtained without the need of forcing a zero value neither for the real 
part nor for the imaginary part of these expressions. 

In Figs. [TU] and ^ we also plotted the planes BR{Kl tt° fi+ h-)-BR{Kl Tr^e+e"), BR{Kl n^e+e-)- 
BR{Kl i:^vv) and BR{Kl ~^ 7r"^+/x~)-BR(ifL tt^vv). It can be seen that the decays Kl n°^+fi~ 
and Kl 7r"e+e~ are subject to weaker modifications than the decay Kl tt^vv. When these plots are compared 
with similar ones obtained in other models, we find clear differences. In this sense, these results can allow to test the 
3-4-1 extension under study here when more accurate measurements of the corresponding decays become available. 
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